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ON THE ACCURACY OF TRIGONOMETRIC INTERPOLATION* 



DUNHAM JACKSON 

Let f (x) be a function of x , with the period 2ir , which takes on known 
values, determined by physical measurement or otherwise, at 2w + 1 points 
equally spaced throughout a period. It is well known that a finite trigono- 
metric sum, of the nth order at most, which coincides with f {x) at the points 
in question, is given by formulae analogous to those which characterize the 
Fourier's development of / ( a; ) , the integrals of the latter development being 
replaced by finite sums involving the known quantities. It has been found 
by Faber that in the case of continuous functions the ordinary sufficient 
conditions, in the way of restrictions on f (x) , for the convergence of the 
Fourier's series are sufficient also to insure the convergence of the interpolating 
function to the value f (x) at all points, as the number n is indefinitely in- 
creased; although in cases not covered by these explicit conditions one develop- 
ment may converge while the other does not.f In the present note it is shown 
that a method by which the author has studied the rapidity of convergence of 
Fourier's series is adapted to the treatment of the corresponding problem in 
interpolation, t and yields similar results. In fact, these results are obtained 
by a simple combination of materials already at hand. It is further pointed 
out how a finite number of observed values may be used to define a formula 
of approximation which converges more rapidly in certain cases than the 
ordinary interpolation formula. 

I. Ordinary Interpolation. 

Let f {x) be a function of x having the period 27r . Let a;i , a;2 , • • • , X2n+i 
be a set of values of x such that 

27r 

Xi+i — Xi = „ I (i = 1, 2, •••, 2ra). 

* Presented to the Society, February 22, 1913. 

t Faber, Mathematische Annalen, vol. 69 (1910), pp. 372-443; see pp. 417-443. 
At the close of the article, Faber considers also certain classes of discontinuous functions. 
The sufficiency of the condition of limited variation had previously been deduced by db la 
Vall^e Poussin, Bulletins de l'Acad6mie royale de Belgique, Classe 
des Sciences, 1908, pp. 319^03. 

t De la Valli^e Poussin (loc. cit., p. 389) obtains the result that if / (a; ) possesses a de- 
rivative having limited variation, the error of the interpolation formula nowhere exceeds a 
constant multiple of 1/w . This theorem is not among those obtained in the present paper. 
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If the coefficients in the trigonometric sum 



[October 



Sn(x) 



+ ai cos X -\- •' • -\- an cos nx + bi sin x + ■ • • -\- bn sin nx 



are determined by the formulae 

O 2ra+l 



ah 



2n 



"TT 2/ (^>) cos Aa;i, 6^ 
~r J- 1=1 



the sum can be expressed in the form* 



O 2»+l 

{h = 0, 1, ■■• ,n), 



Zn+l 



sin [ ( ra + J ) ( a;,- — a; ) ] 



''"^^^ 2n+l£?''^^'' sin 1 (a;.- -a;) 

by means of a famihar identity, and from this expression it appears that 
Sn (xj) = / ( Xj ) at each of the 2ra + 1 points Xj , all the terms but one under 
the sign of summation being zero. It is the only trigonometric sum of order 
n or lower which has this property, since two such sums which are equal for 
more than 2ra distinct values (mod 2t ) of the argument must be identical. 

It was proved by FASERf that if / ( a; ) never exceeds 1 in absolute value, 
then the absolute value of S„ (x) can not exceed a certain constant multiple 
of log n . The proof, which may be given here for the sake of completeness, 
is as follows. With the hypothesis made concerning f (x) , 



Sn{x) 



< 



2n + 



2TC + 1 

J- i=l 



Set ^ {Xi — x) 



Ui] then 

\Sn{x)\ 



^ 1_ 

< 2n-\- 



sin[ (n + h) (xj — x)] 
sin ^ {xi— x) 

'"+' sin(2n+ !)«.■ 



1^ 



sin Ui 



In view of the periodicity of the quotient, it may be assumed without loss of 
generality that all the variables m, lie in the interval from — ^tt to | tt . Let 
t'o be the smallest of the numbers | Ui \ (or one of the two smallest, if two are 
equal), Vi the second, and so on. Then we have, on rearranging the sum. 



I S„ {x) 
where, as is readily seen. 



= ^ sin(2n+ \) v, 
j=o (2?i + 1) sint)i 



ITT 



Now 



= ^ (i+l)7r 

2(2n+ 1) < ""'■< 2(2ra+ !)■ 



* Throughout this paper, an expression which, regarded as a quotient, becomes meaningless 
for certain values of the variable by taking on the form 0/0 , will be understood to be defined 
for those values so as to be a continuous function. 

t Loc. cit., pp. 419-420. 
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(1) 



sin ( 2n + 1 ) »o I | sin ( 2n + 1 ) «o 



Vo 



^l, 



( 2n + 1) sin Vq \ ! (2n + 1 ) »o sin Vq 

since ( sin » ) / » is a monotonically decreasing function in the interval ( , -tt ) 
which includes both the points Vo and ( 2n + 1 ) »o • Similarly, the absolute 
value of the term in Vi does not exceed 1. For i^ 2 , 



sin ( 2n + 1 ) Vi 



( 2n + 1 ) sin vi 



< 



(2n+ l)sin»,' 



s'nce sin » > 2» / x for < » < Itt , we may write 

(2) (2w+ 1) sinsi ^ (2?i+ 1) sin 

and hence 



IX 



2(2w+ 1) 



> i, 



^ sin(2w+ i)vi 
t=2 ( 2?i + 1 ) sin Vi 



2n 1 

^ Sj< log2n. 

i=2 '' 



Consequently 

I S„ (a;) I < 2 + log2 + logn< 51ogn (nS2). 

It would be easy to replace the constant 5 by a smaller value, but we shall not 
make any effort to do so here. 

From the inequality just obtained, it follows immediately that if | / ( a; ) | < e, 
then I iS„ ( a; ) I < 5€ log n , and \f {x) — S„ (x) \ < €-{• 5e log n <_ 7e log n , 
nS 2. 

Now let the restriction on f {x) be replaced by the broader hypothesis 
that there exists a trigonometric sum Tn {x) , of the nth order at most, such 
that for all values of x we have the inequality 

\f{x)- Tn{x)\ ^ €. 

The interpolating function Sn{x) that corresponds to / ( a; ) may be obtained 
by adding those for Tn{x) and f or / ( a; ) — Tn {x) . The former is identical 
with Tn {x) , being equal to it at 2n + 1 points, and so an error is committed 
only in approximating to the difference f (x) — Tn (x) , which does not 
exceed e. It follows that the error can never be greater than ?€ log n if 
n ^ 2 , and we have, to sum up,* 

Theorem I. 7// ( a; ) can be approximately represented by a finite trigonometric 
sum of the nth order at most, n ^ 2, with an error never exceeding e , then f (x) is 
represented by its interpolating function Sn{x) with an error not exceeding It log n . 

This can be immediately combined with theorems of the author concerning 
the accuracy of approximation by trigonometric sumsf to give a variety of 
more explicit results. To mention only two: 



* Cf. Lebbsgue, Sur les intigrales singulieres, Annales de la Faculty de Tou- 
louse, ser. 3, vol. 1 (1910), pp. 25-117, for the corresponding result in the case of Fourier's 
series (see pp. 116-117). 

t These Transactions, vol. 13 (1912), pp. 491-515. This article will be referred 
to hereafter as A. Also, these Transactions, pp. 343-364 of the present volume. 
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If lims=o [f (3; + 8) — f (x)] log I 5 I = 0, uniformly for all values of x, 
then Sn (x) uniformly approaches f (x) as its limit when n is indefinitely 
increased.* 

If /(a;) everywhere satisfies the Lipschitz condition f 

\f(x")-f{x')\^\\x"-x'\, 
then, for all values of x , and f or n > 2 , 

The preceding work also makes it easy to estimate the possible effect of 
errors I in the determination of the quantities /( aij ) . Let Sn (x) be the 
interpolating function defined by the true values f (xi) , and Sn (x) the func- 
tion defined by a set of observed values J (xi) . The difference between 
these two functions is the corresponding function defined by the differences 
/ ( Xi )—]{Xi), and, if the greatest of the latter does not exceed e , the maximum 
of \Sn{x) — Sn{x)\ can not exceed be log n . It would be possible to replace 
the coefficient 5 by a smaller one, as has been said before, but it is interesting 
to remark that it would not be possible to dispense with the factor log n, 
or even to replace it by one which becomes infinite less rapidly with n , since, 
as Faber shows, the maximum of | Sn ( a; ) | for functions f (x) not exceeding 
1 in absolute value actually is of the order§ of log n. 

II. A Modified Formula. 
For comparison with Sn (x) , let us define a function S^ ( a; ) as follows : 

(3) ^-M = H4fM[^^i^] 

where 

i_ = y [ sm^mjxi- x) Y 
Hm i=i\_m sin ^ {xi — x) \ ' 

The points Xi are assumed to be such that|| 

* This was deduced by Faber (loc. cit., p. 422); the particular theorem used here concerning 
approximation by trigonometric sums was proved by Lebesgue, loc. cit., pp. 115-116. 
t See A, Theorems I, VI. 

I The fact that, in the case of polynomial interpolation with equidistant ordinates, errors 
of observation may produce a disproportionately large effect, was pointed out by de la 
Vallee Poussin, loc. cit., pp. 321-322; see also p. 346 of the same article. 

§ Faber, loc. cit., p. 424. If / ( Xi ) = ( — 1 )*"' , and ii — x = ir/C 2ra + 1 ) , then 

/ ( Xi ) sin [ ( n + i) ( Xi - X ) ] = + 1 
for all values of i , and it is readily seen that S„ ( x ) is of the order of log n . 

II It would be possible, without essential change, to use instead an odd number, 2to — 1 , 
of equally spaced points. 
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It 
Xi+i — Xi = — ( i = 1 , 2 , • • • , 2m - 1 ) . 

m 

Let us first notice that the apparent dependence of Hm on the variable x is 
illusory; Hm really depends only on m. This may be seen as follows. The 
fourth power under the sign of summation can be expressed in the form* 

AmO + Ami cos {Xi — x) + Ami COS 2 (Xi — x) + • • • 

+ Am, 2m-2 cos [ ( 2m — 2 ) ( a;, — a; ) ] 



2m— 2 

E 

A=0 



= 2 Amh ( COS hxi COS hx + sin hxi sin hx ) , 



the coefficients Amh being constants. It is a well-known property of the tri- 
gonometric functions that if 2/1 . • • • . Vq are q points disposed at successive 
intervals of 2Tr / q, and p is not divisible by q, thenf 



1 9 



23 cos pyt — 2 sin pyt = . 



<=i *=i 



For the application that we wish to make, q is 2m , and p is one of the numbers 
, 1 , • • • , 2m — 2 , and so the sum defining 1 / Hm reduces to 2m^mo . 

From this fact and the expression used above for the fourth power involved, 
it appears that Sm ( a; ) is a trigonometric sum in a; , of order not higher than 
2(m— 1). It is a formula of interpolation, in the sense that it is completely 
determined by a finite number of values of the function represented, though 
it does not, like (S„ ( a; ) , become equal to / ( a; ) at the points Xi , in spite of the 
fact that it is of higher order than /S„ (a;) in proportion to the number of 
points used. We shall proceed to establish the property which constitutes 
its chief claim to attention: If /(a;) satisfies a Lipschitz condition, it is re- 
presented by Sm (x) with an error which approaches zero at least as fast as 
1 / n in order of magnitude, if we denote by n the maximum order 2m — 2 
of the trigonometric sum. 

Suppose that f (x) satisfies everywhere the Lipschitz condition 

|/(a;")-/(a;')|^X|a;"-a:'|. 



* See A, p. 493. The absence of sine-terms follows from the fact that the expression is 
an even function oi Xi — x . 

t Cf . e. g. B6cHER, Introduction to the theory of Fourier's series, Annals of Mathe- 
matics, ser. 2, vol. 7 (1906), pp. 81-162; p. 135. For an analytic proof, a starting point 
is the fact that when p = 1 and one of the yt's is zero, S cos yt and 2 sin j/i are the real and 
imaginary parts respectively of the sum of the gth roots of unity, which is zero if g > 1 . 
Then the cases may be considered successively that p is any number relatively prime to q , 
and that p and q have a greatest common divisor greater than 1 but less than q, still with the 
assumption that one of the yt's is zero. The last generaUzation may be effected by letting 
2/t = 2/1 + yi' , where 2/1 is arbitrary. 
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We may write 

or, setting ^ (Xi — x) = ut, 



\i:^(x)-fix)\ = 

(5) 



ff.|[/(. + 2„.)-/Wl[|^]*| 

7^1' 'Lmsinwij 



Because of the periodicity of the functions involved in the first line of this 
last relation, we may assume without loss of generality that all the variables m,- 
there, and so also in the following line, are in the interval from — ^tt to + 5^ . 
The work of determining approximately the magnitude of Hm has essentially 
been performed in the article A. For, if we recall the notation Jm of that 
article,* we have 

r'^ r sin mw "I ^ , C"'^ , , mir , ir 1 

Jo L ™ sm M J Jo 2 i Hm 

As it was shown that ^ir ^ Jm > iir , for all values of m , it follows that 

1 ^ 77 ^ s. 

2 <. -"m '^ 4 • 

In the sum which still remains to be evaluated, let the smallest of the quan- 
tities I Ui I (or, if this is indeterminate, one of the two smallest) be set equal 
to vq , the next to vi , and so on. Then 

^ , , r sin mui V '^' f sin mvi I* 
7^5 ' ' L m sm Mi J S L ^ sm »»• J 
where 

iw — _ ( i + 1 ) TT 

Adapting the relation (1) to the present case, we see that 

[sin mvo "I * == , f sin mvi "1 * ^ ^ 

msin»oJ ' Lmsini^iJ 

For i ^ 2 , by reasoning similar to that employed in connection with (2), 



and 



Hence 



_ .117 — 

m sm Vi ■> m sm -. — > 



r sin mvj "1 ^ = /2Y 
\_ni sin Vij \i / ' 



' A, pp. 503-507. 
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Now »o + ill = T / ( 2m) , as is seen by reference to the definition of the symbols 
Vi-, and 



du\ 



~ m \2'''3/ 3m * 
Consequently 

'y^' r sin mr» T a- / 1 10\ 237r 
^ * L wi sin 1),- J m \ 2 3 / 6m ' 

and so, by (5) and the inequality for Hm , 

If we denote by n , as was suggested before, the maximum order of 2^ ( a; ) , 
which is less than 2m, we have finally: 

Theorem II. If the function f (x) , of period 2ir , everywhere satisfies the 
Lipschitz condition | / ( a;" ) — / ( a;' ) | ^ X | a;" — a;' | , it is approximately 
represented by the finite trigonometric sum 2^ ( a; ) , of order not higher than n , 
with an error not exceeding 37X / n . 

The coeflBcient 37 is not particularly significant; no attempt has been made 
to reduce it to as small a value as possible. The main point is that we have 
found here an upper limit for the error which is of higher order in its approach 
to zero than that which we had for the ordinary interpolation formula. We 
have not proved, to be sure, that the order of the limit of error which we found 
in the earlier case is the best that can be obtained, but it is not difficult to 
show, by the method applied by Lebesgue* to the corresponding problem in 
the case of Fourier's series, that the error of the approximating function 
Sn (x) for a function f (x) satisfying a Lipschitz condition can actually 
exceed a constant multiple of ( log n) j n for infinitely many values of n . 
The expression 2^ ( a; ) accordingly is really superior to the other in the order 
of its approximation to certain functions of the class considered. 

It may be observed further that errors in the quantities f (xi) , none of which 
exceeds e in absolute value, can not affect 2m ( a; ) at any point by more than 
€ , as compared with a quantity of the order of e log n in the case of tS„ (x) , 
and that 2^ ( a; ) converges uniformly to the value / ( a; ) , as m is indefinitely 

♦Lebesgue, Bulletin de la Soc. Math, de France, vol. 38 (1910), pp. 184- 
210; pp. 203-206. 
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increased, for every continuous function f (x) . The former assertion is 
coextensive with the statement that | S^ ( a; ) | is never greater than 1 if 
I / ( a; ) I never exceeds 1, a fact which is obvious from the definition of Sm ( a; ) . 

To prove the convergence property suppose that e is an arbitrary positive 
quantity. Let 5 be a positive number such that \f(x") — f (x') \ < ^e 
whenever \ x" — x' \ < 5; the existence of 5 is a consequence of the hypo- 
thesis that f (x) is everywhere continuous. Let x have any particular 
value, and divide the sum on the right hand side of (4) into two parts, 
one containing the terms* for which {xi—xl'^S, and the other, the 
remaining terms. The first part does not exceed ^e in absolute value, for 
any value of m , in consequence of the definition of 5 . In the second part 
of the sum, \fixi) — f (x) \ does not exceed twice the maximum of \f (x) \ , 
and sin ^ (xi — x) ^ sin J 5 ; consequently each term of the second part is 
inferior in absolute value to 1 / m* multiplied by a quantity independent of 
m and x, and the whole second part does not exceed a constant multiple of 
1 / m' . As soon as m is greater than a suitable quantity, independent of x , 
the second part too becomes smaller than |e , and the uniform convergence is 
established. 

The properties of S^ (a;) last developed, while interesting in themselves, 
are less characteristic than the one on which the principal stress has been laid; 
they are shared, for example, by the simpler expression which is obtained if 
the fourth power in the definition of S,„ ( a; ) is replaced by the corresponding 
square. 

To apply the formula Sm (a;) to the problem of polynomial interpolation, 
the given ordinates, to be sure, being unequally spaced, suppose that g (y) 
is a function of y defined in the interval — 1 < ?/ < 1 • Suppose that the 
values of g (y) are known, either at the m points 



(2i- l)7r 
yi = cos 

or at the m + 1 points 



2/i = cos ^ (i = 1, 2, •••,to), 

2m 



2/» = cos — ( i = , 1 , 2 , • • • , TO ) . 

m 

If we set y = cos x, g (y) = g ( cos x) = f (x) , we have the values of f (x) 
given for a set of values of x symmetrically situated with respect to the origin, 
and equally spaced at intervals of tt / m . Let those which lie in the interval 
— ir < a; ^ IT be used as the points Xi for defining a function Sm (a;) cor- 
responding to our present function f (x) . If we replace a; by — a; in this 
expression, and at the same time replace Xi by — a;i , which amounts only to 



* We suppose the numbers Xi replaced, if necessary, by values congruent to them (mod 2ir), 
so that none of the differences \xi — x\ is greater than x . 
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a rearrangement of the terms in the sum (3), because of the symmetrical 

arrangement of the points Xi , we see that 2^ ( a; ) remains unchanged. Being 

an even function, it is a polynomial in cos x , of degree 2 ( m — 1 ) at most. 

We have thus a means of approximating to gr ( ?/ ) by a polynomial in y . The 

properties of the trigonometric interpolation formula developed above can be 

immediately interpreted with reference to the present case.* 

Harvabd Univebsity, 
Cambridge, Mass. 

* Cf . A, p. 495. 



